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ABSTRACT 


Ponty two ocean bottom gravity stations in southern 
Monterey Bay were occupied in the summer and fall of 1972 from 
the R/V ACANIA. A land gravity survey of ten stations about 
the perimeter of the Bay was conducted in the se οὐ ΠΟ. 
Gravimeters employed were LaCoste and Romberg Models H6G and 
G-17B, respectively. 

Conventional steps in data reduction are discussed, and a 
terrain correction theory unique to ocean bottom gravimetry 
is presented. The complete Bouguer anomaly (CBA) field for 
bottom and shoreline surveys is included. 

The geological interpretation of the gravity data is 
discussed briefly.  Sub-bottom structure of southern Monterey 
Bay as determined by seismic reflection is verified by the 
CBA field, and a calculated density contrast between the 
basement granodiorite and overlying sedimentary strata is 
Bound to be realistic. The data supports the existence of a 


fault oriented beneath the Monterey Submarine Canyon. 
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IN RODUGTION 


/ 


Geological oceanographers have three principal means of 
investigating structures below the ocean floor: seismic, 
magnetic, and gravimetric surveys, the latter of which is the 
Sue? of this writing. The relatively small variations in 
the acceleration of gravity over the surface of the earth are 
largely due to such obvious factors as position, elevation, 
crustal density, and local terrain; once these effects which 
may be unrelated to substructure irregularities are removed, 
however, the investigator is left with a meaningful residual 
in the form of the complete Bouguer anomaly (CBA). This mea- 
sure of anomalous attractive force per unit mass then consti- 
tutes a basis for the interpretation of the upper few hundreds 
or thousands of feet of the earth's heterogeneous crust. 

This study is based upon measurements made at 10 land 
stations and 82 underwater stations which were occupied in the 
spring and summer of 1972. The bottom gravity data collected 
was obtained aboard the Naval Postgraduate School's (NPS) re- 


search vessel R/V ACANIA (Fig. 1). 


AS ETE Y PURPOSE 

The main purpose for conducting this survey was to learn 
more about the structural geology of southern Monterey Bay by 
obtaining gravity data in areas too shallow for larger vessels 
capable of surface gravimetry. Such vessels must stabilize 


their course and speed for several minutes after each turn 
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ene useful data can be obtained; this usually places pro- 
hibitive restrictions on how close to shore the vessel may 
work, depending upon local bathymetry. But such geographical 
Actions do not apply to bottom gravimetry, since the 
research vessel can slowly decelerate to even very shallow 
stations and stop for meter lowering and measurement. The 

R/V ACANIA was able to survey within 100 yards of the beach in 
water less than 50 ft deep. Additionally, bottom gravimetry 

is generally more accurate than sea surface gravimetry due to 
the stability of the measuring platform, the sea floor.  An- 
other advantage of bottom gravimetry is its closer proximity 
meme density contrasts that produce a gravity anomaly; sur- 
face ship meaSurements must contend with the attenuating effect 
Of vertical distance from the meter to the Nc Uc (Beyer, 
et al., 1966). This means that the deeper the station, the 
greater the accuracy of bottom gravimetry over surface gravi- 
metry (assuming equal meter precision). 

After a brief area orlentation, this report will describe 
some basic theory of gravity measurements, specific details on 
nent used, and various procedures employed during the 
survey itself. This will be followed by a section explaining 
each step in the reduction of the raw data, and an analysis of 
the results in terms of geological significance will conclude 


τας report. 


ΠΕΕΑ DESCRIPTION AND GEOLOGICAL SETTING 
Monterey Bay is located about 70 miles south-southeast of 


Sr aneilsco, California. It is characterized by gently 
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sloping sandy bottoms in the north and south, divided by one 
of the largest submarine canyons in the world. The area 
encompassed by the author's survey includes that portion of 
southern Monterey Bay within the 50-fm depth contour from 
Pinos to Moss Landing (Fig. 2). 

MEN the exception oí the rocky bottom in the extreme 
southwestern portion of the survey area and the predominance 
of green muds near the Monterey Submarine Canyon, a sandy 
ocean floor exists in southern Monterey Bay. Cretaceous 
Santa Lucia granodiorite, readily visible at Pt. Pinos and 
Eu crtuting much of the adjacent sea floor, is buried to in- 
creasing depths beneath less dense sediments and sedimentary 
rocks as one proceeds toward the Salinas River mouth-Moss 
Landing area (Fig. 3). Much of the overlying material to the 
north and east is the Monterey Formation, a Miocene rock 
consisting largely of silicious mudstone, diatomite, and marine 
Shale, interbedded with beds of opaline chert. Above the 
Monterey Formation are thin layers of the Pliocene Paso-Robles 
Formation and Pleistocene Aromas red sands.  Thinner layers of 
A kaic deposits from the Salinas River have been laid down 
above these formations during the Holocene period. Additionally, 
the Pliocene Purisima Formation, abundant in northern Monterey 
Bay, has recently been found along the south wall of the 
Monterey Submarine Canyon (H. G. Greene, personal communication, 


973). 
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Figure 3. Geological Map Of the Survey Area 
(after Greene, 1970). 
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eee PREVIOUS WORK 

Past geological investigations of Monterey Bay have 
mainly concerned themselves with the origin and present 
structure of the Monterey Submarine Canyon, as seen in the 
studies by Shepard (1948), Martin (1964), and Martin and 
Emery (1967). Geological research in the shallow water 
regions of the Bay has been conducted, but most of these 
efforts investigated specific geological properties, such as 
bottom sediment type and distribution. The first extensive 
and detailed report of the geology of southern Monterey Bay 
was recently presented by Greene (1970). This work reviewed 
the results of previous bathymetric data, rock dredges, and 
grab and core samples, as well as the author's own seismic 
reflection survey results, in order to make inferences con- 
cerning the salt water intrusion problem in the aquifers of 
the lower Salinas Valley. 

While gravity data has heretofore been totally absent in 
E-nallow off-shore regions, various investigators have 
Surveyed adjacent land areas. Fairborn (1963), Sieck (1964), 
and Ivey (1969) have surveyed the local Monterey-Salinas-Ft. 
Ord region, and Bishop and Chapman (1967) compiled land 
gravity data from many sources in preparing their Bouguer 


gravity map of the area. 
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Ην GRAVITY MEASUREMENTS 


D. INSTRUMENTS 

When the variation of a quantity is extremely small com- 
pared to the value of the quantity itself, measurement of that 
variation requires some method of exaggeration in order to 
achieve desired sensitivity. Such is the case at hand: 
gravity values at the equator and the poles are approximately 
978 and 983 gal, respectively (owing mainly to the poleward 
decreases of the local radius of the planet and the centrif- 
ugal force due to its rotation); hence, the entire range of 
possible variation relative to the mean value is about 0.53. 
Additionally, this small total range of values is about a 
hundred times greater than the range observed in Monterey Bay, 
SO the requirement for precise instrumentation is clear. Most 
of the numerous different gravimeter designs that have become 
operational over the years illustrate this idea. Depending 
upon operating principle, they fall under one of two general 
categories: stable ana unstable systems. 

[ieee stable Gravimeters 

Stable gravimeters use a spring-balance system wherein 

the attractive force on a known mass suspended from the end of 
a spring is linearly proportional to the resultant elongation. 
In the case of the Hartley gravimeter, the minute elongation 
is made more measurable by passing a small, horizontally- 


hinged beam through the spring-mass system and attaching a 





mirror close to the unhinged end, vhere beam displacement is 
greatest (Fig. 4). Light from a simple remote source is re- 
flected onto a graduated scale for the observation of relative 
differences. A micrometer screw attached to an auxiliary 
Spring is located near the hinged end of the beam so that the 
light beam may be brought Lack to some fixed reference position. 
The amount the micrometer screw is turned to accomplish this 
zeroing process is thus a measure of the difference between 
the local gravity value and the value associated with the re- 
ference position. The system 1S termed stable because the only 
ea acting to counter gravity is the spring restoring force. 
2. Unstable Gravimeters 

For deflections of the mass in the unstable (or astatic- 
balance) gravimeter, the force of gravity is supplemented in 
Mmeseasing the spring elongation's restoring force by a third 
force which acts in the same sense as that of gravity itself. 
The Thyssen gravimeter iS a typical example (Fig. 5). The 
clockwise moment of an auxiliary mass above the pivot point 
results in additional beam displacement which can be kept 
ENICar!y proportional to gravity. Thus, greater sensitivity 
is the chief advantage of the unstable gravimeter; for a given 
change in gravity, the unstable device will undergo greater 


beam displacement than that of a comparable stable gravimeter. 


B. OCEANIC VERSUS LAND GRAVIMETRY 
In the case of land gravity meters, systems like those 


Mescribed above are enclosed in an instrument roughly the size 
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Figure 4. Schematic of a Stable Gravimeter: 
the Hartley Model (after Dobrin, 1960). 
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Figure 5. Schematic of an Unstable Gravimeter: 
the Thyssen Model (after Dobrin, 1960). 
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of a kitchen toaster, and are therefore as mobile as the 
operator. The mass is manually clamped in place after each 
reading and the only electrical requirement is a battery to 
pusude power for the light source, scale illumination, and 
temperature regulation (necessitated by the temperature- 
dependence of the expansion coefficients of the springs). 

Dich regard to oceanic gravimetry, however, the situation 
is vastly more complicated. Three platform options are: 
Surface ship, submarine, and bottom measurements; each has 
its own drawback. Surface ship gravimetry requires intricate 
sensors to measure the effects of ship-motion accelerations. 
Submarines designed for See eee ee research are still quite 
rare; military submarines have severe space limitations and 
the more pressing function of national defense to attend to. 
Bottom gravimetry (conducted with operators accompanying the 
es the bottom in diving bells prior to the 195015) in- 
volves the expensive business of remote control. In addition 
to these monetary problems, all types of oceanic gravimetry 
must contend with the formidable problem of accurate navigation, 
a difficulty not common to most land surveys. 

Other difficulties are inherent and peculiar to bottom 
gravimetry. The gravimeter itself must be enclosed within 
some type of waterproof chamber strong enough to withstand 
pressures of many atmospheres. Also, a large number of 
separate electrical conductors, normally within the insulation 
and armor of an oceanographic cable, must be employed to effect 


the numerous recordings and operating modes. This requires a 


zu 





EpSsemalized winch with a multitude of slip-rings. Rocky or 
Steep bottoms can incline the instrument to an extent beyond 
the range of the leveling gimbals. But perhaps the most 
challenging problem of all is the requirement that the re- 
Search vessel maintain its position above the stationary 
EI despite drift induced by wind and current. If this 
requirement is not met, and the winch operator has not lowered 
sufficient excess cable, the meter may be dragged and jarred 
while in the unclamped mode, thereby damaging the meter's 


delicate spring suspension system. 


22 








tite. EQUTPMENT 


A. GRAVIMETERS 

Of the two general types of meters discussed in the pre- 
vious section, the LaCoste and Romberg models fall in the un- 
Sable category. Figure 6 1S a drawing of the main components. 
When the mass is released from its protective clamps for the 
moose of taking a reading, its motion to an equilibrium 
position is impeded by very heavy air damping; it therefore 
takes a finite time interval for the mass, support beam, and 
light beam to come to rest. The speed of the motion of the 
light beam is a function of the position of the measuring screw 
and 1S indicated by a beam position galvanometer which may be 
Esnnected to a strip chart recorder. The slope of this record 
of time versus speed is a measure of how close the measuring 
ED (SS to the equilibrium position, which, in turn, is a 
measure of gravity. 

The zero-length spring in Figure 6 was first introduced by 
LaCoste in 1934 for use in an innovative Tone period. vertical 
seismograph. Its unique method of manufacture renders elon- 
gations or compressions proportional to increases or decreases 
Io ravity itself, resulting in increased sensitivity. 

The calibration of most other types of gravimeters depends 
on weak auxiliary springs which are connected directly to the 
measuring screws (Figs. 4 and 5). This is avoided in LaCoste 
and Romberg gravimeters by employment of a lever system 


connected to the measuring screw. The light beam is thereby 
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Figure 6. Simplified Diagram of the LaCoste 
- and Romberg Gravimeter (after LaCoste, 
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上 ei to the equilibrium position by actually changing the 
location of the Zero length spring's- upper support. 

l. Lacoste and Romberg Model G-17B Geodetic Land Gravimeter 

The LaCoste and Romberg model G-17B geodetic land 

gravimeter used for the shoreline survey is an instrument 
housing the system just described in a 5-3/4 inch x 6-1/2 
1409 6-1/2 inch case weighing about 5 Ib. This model was 
mE ntroduced in 1956, and operates in the field on battery 
power over a 7000 mgal range. Instrument drift is less than 
0.5 mgal/month due to a thermistor-transistor heater control 
system. Sealed against atmospheric pressure changes, the G-17B 
gravimeter is nulled by adjusting the measuring screw while 
observing the light beam positon through a microscope eyepiece. 

2. LaCoste and Romberg Model H6G Underwater Gravimeter 

The LaCoste and Romberg Model H6G underwater gravimeter 

used for the bottom survey is a complete system which permits 
mStrument like the G-17B to be used on the ocean floor. 
eemply stated, 1t is a remote-controlled land gravimeter in- 
side a shell of two thick aluminum hemispheres supported by a 
angular Base (FIS and aA brief description of the 


features of this model is presented in Table I. 


BSEAUXTETARY EQUIPMENT 

The bottom gravimeter and all of the associated equipment 
were loaned to the NPS Oceanography Department by the U. S. 
Naval Oceanographic Office in Suitland, Md. Engineered and 
assembled by the Geophysical Equipment Manufacturing Co. of 


Houston, Tex., the auxiliary equipment's operation may be 
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Figure 7. sodel H6G Gravimeter Ready for Use. 
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Table I. Features of the LaCoste and Romberg 
Model H6G Underwater Gravimeter 


Base Size 


Dinei egui lateral triangle 
(with legs) B d 


Height 29 inches 


Weight 350 lb, including enough weight 
to insure acceptable sink rate 


Optimum Accuracy πο ue に 0l 
Realistic Accuracy +0.10 mgal 


Routine Maximum 















Operating Depth 99. ene 
Maximum Operating 

Depth Tested ο 
Maximum Bottom 15? 


Slope Permissible 


Pressure sensor (0.5% accuracy) 
connected to casing port 





Depth Indication 








Optically-controlled servo motors 


Leveling System driving rack and pinion assemblies 


Open circuit, 1/4-inch gap, inside 
bottom of lower hemisphere 


Flood Indication 











Summarized as follows: a gasoline engine provides mechanical 
energy to a hydraulic pump transforming the rotary motion into 
hydraulic pressure (Fig. 9). This hydraulic pressure is used 
to operate both the winch and A-frame using two, two-way con- 
trol valves. One valve controls a hydraulic motor which con- 
verts the fluid pressure back into rotary motion to drive a 
heavy-duty chain link connected to the winch reel's shaft. 
The A-frame is hinged at deck level and actuated from above by 
two hydraulic pistons controlled by the other valve. 

With the gravimeter connected electrically to the termina- 
INE USt one end of the cable, the other end is routed inside 
the hollow shaft of the winch and out the side to a slip-ring 
assembly. The outer side of the slip-ring assembly is wired 
mena cannon plug connection on the control box, which in turn 
15 connected through a rectifier and an isolation transformer 
Pone of the ship's 115 vac outlets, the basic power source 


for the entire electrical system. 
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Figure 9. Schematic Diagram of the Auxiliary Equipment. 
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I, SE EROCEDURE 


This survey was part of a larger effort to investigate 
off-shore gravity values over a more extensive area: two co- 
workers (also NPS students) recorded data in northern Monterey 
Bay and Carmel Bay. Since bottom gravity studies are by no 
means a one-man operation, they were assisted by the author 


and provided similar assistance to him during his survey. 


A. GRAVIMETER CALIBRATION 

Prior to beginning the shoreline survey on 26 April 1972, 
the land meter was calibrated by the author and co-workers 
under supervision of U. S. Geological Survey (USGS) personnel. 
Readings were taken at several benchmarks from USGS head- 
quarters in Menlo Park (USGS 1 JD) to Skeggs Point Scenic View 
(meso. Coast and Geodetic Survey Station B-388), a frequently 
used calibration route spanning a range of 137.20 mgal (Chapman, 
1966). Subsequent reduction of counter values to gravity values 
resulted in a difference of 137.15 mgal over the entire range. 
Smaller differences were found at the intermediate benchmarks. 

The underwater gravimeter was initially calibrated by the 
manufacturer prior to shipment to Monterey in the spring of 
1972. After the survey was completed, another calibration 
check of the meter was in order. Early in 1973 the author 
Occupied CA-259 (Woollard Airport Base WA-84) at the Monterey 
County Airport and a benchmark located at the base of the steel 


tower on the Monterey Coast Guard pier, designated WH-29 by 
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Woollard and Rose (1963). The meter readings resulted in a 
gravity difference less than 0.5 mgal greater S NETS 20955 
mgal established value (Chapman, 1966). The difference was 
probably due to the recent construction of a second terminal 
IEuEdIng next to CA-259. Since the airport value is lower 
than the pier value, and since the effect of the new building 
is an additional upward attraction, the difference between the 
two reference stations has probably increased due to the new 
Gemetruction. In view of this consideration, as well as the 
E x:that the published location of the CA-259 benchmark could 
nube found exactly (due to the construction), it is estimated 
that the observed gravity difference between the airport and 


pier stations was within 0.1 mgal of the true value. 


ιν... SHORELINE SURVEY 

A preliminary gravity survey of the perimeter of southern 
Monterey Bay was undertaken in late April, 1972, with the land 
gravimeter discussed above. The information from this shore- 
EE vey of 10 stations was later used to verify the near- 
by CBA values and isoline trends of the underwater survey, 
Since land gravimetry is an order of magnitude more accurate 
than oceanic gravimetry. (While land gravimeters are read to 
the nearest 0.01 mgal and subsequent corrections' accuracies 
can lead to CBA values of 0.1 mgal accuracy, remote-controlled 
bottom gravimeters, read to the nearest 0.1 mgal, can result 
Ma. CBA accuracy of at best 1 mgal.) In addition to verifying 


the bottom survey results, these land stations also served to 
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tie in the combined survey with previously surveyed areas in- 
land. All 10 land stations were located (vertically) within 

a few feet of sea level so that the small magnitude of the 
Subsequent gravity corrections would yield better CBA 
accuracies. The stations were fairly evenly spaced from Pt. 
Joe to Moss Landing (Fig. 10). Visual estimates of station 
level relative to the average existing sea surface elevation 
were recorded at each station, so that station elevations rel- 
ative to mean sea level could be calculated later from tidal 
information. Exact locations and depths/elevations of both 


underwater and shoreline stations are included in Appendix A. 


Seok BOARD INSTALLATION AND SEA TRIALS 

When the auxiliary equipment arrived in Monterey in late 
June, 1972, preliminary inspection indicated that it had been 
some time since the system was used. A re-termination of the 
cable and an engine overhaul were required before the survey 
"md begin. 

An initial concern was the question of location of the 
equipment aboard the ship. It was decided that the after 
section of the upper level would be most suitable in view of 
the physical requirements associated with meter lowering and 
the intended length of stay of the equipment aboard. The A- 
frame, supporting plate, and gravimeter were stationed on the 
Starboard side so that the data recorder in the dry lab directly 
below would be able to observe cable speed and direction of 


Doelon (Fig. 11). Also, the crane's location in close proximity 
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on the starboard side seemed compatible with the plan to NU. 
meeter to and from the main-deck for each station. The 
winch was placed abeam the A-frame area but on the port side, 
ο οεεεπρε to spread out the 4 ton weight of the equipment. 
The engine and hydraulics stand were located on the ship's 
centerline and slightly forward of the other components. 

A brief sea trial period was arranged prior to the start 
Of the actual survey in order to acquire a familiarity with 
Semercions to be encountered. A trial station was occupied 
west of Seaside, California, in about 100 ft of water. SCUBA 
divers accompanied the meter to the bottom to observe the 
Sink rate, take some underwater photographs, and determine if 
excessive cable would foul the feet of the gravimeter.  How- 
ever, a strong wind and too little excess cable resulted in 
dragging the meter across the sandy bottom while the meter 
waS in the unclamped mode. It was decided to return to the 
ACANIA's mooring to check for possible damage to the mass 
ENSpensron system. Fortunately the counter reading reproduced 
the value obtained just prior to departure for the station, a 
sign cT no damage had been sustained. 

This preliminary venture into the Bay showed that the 
Stations would have to be occupied with the ship's bow into 
the wind so that the major wind-induced position corrections 


could be accomplished with ahead or astern engine commands. 


Since the ACANIA has no bow-thruster, the bow could be kept 
into the wind with alternate use of one or both of the engines 


and rudder. Also, a feel for the amount of cable payed out as 
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a function of cable angle and tension was obtained, since the 
divers reported that cable fouling about the base of the 
meter was a remote possibility (owing mainly to the cable's 


large diameter and resultant large radius of curvature). 


EM SURVEY OPERATIONS 

τος το the onset of the survey, selection of station 
locations was governed mainly by the requirement for a fairly 
evenly-spaced grid from which an accurate scalar field of the 
CBA could be extracted. In view of the size of the survey 
area and the expected ACANIA availability periods, it was felt 
that a minimum of 75 stations would be needed to satisfy this 
requirement. A chart showing the spacing of about 100 intended 
stations similar to Figure 10 was provided to the bridge for 
course information. Whether the actual stations occupied coin- 
cided with these intended positions was not of great concern; 
the important thing was to determine the position of each 
station as accurately as possible during and immediately after 
mee Cupation in order to adjust the planned location of the 
next station for a possible gap in coverage. Stations were 
numbered in order of occupation (except for those on land). 
Upon completion of the survey, the 92 stations occupied covered 
an area of 50 sq n mile, representing a station density of 
almost 2 per sq n mile (Fig. 10). 

I Navigation 

With the exception of three stations occupied during 


morning fog (and consequently requiring radar positioning), 
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all of the bottom station positions were determined by visual 
lines of bearing from prominent landmarks during hours of both 
daylignt and darkness. Three or four such bearings were taken 
at each station, and these were often supplemented by radar 
distances from coastline features. Standard navigational 
procedures were used to estimate the true location of the 
ship's position in the small area where lines and arcs of 
position converged. Most of the three-line position triangles 
were on the order of 0.1 n mile or less on a side, and a good 
portion of the plots were point fixes.  Bearing-taking was 
begun as soon as the ACANIA became dead-in-the-water (before 
the meter was lowered to the bottom) and was always completed 
well before the meter was raised after the measurement. When 
it appeared the vessel was drifting, a second set of bearings 
was taken. Differences in meter position and ship's position 
at each station were thus small enough to safely neglect 
consideration. In view of these factors, it is stated with 
een degree of reliability that station position errors at 
Sea were at worst 0.l n mile. If such an error was in a 
north-south direction the value of theoretical gravity and 
Seema! CBA would differ by only 0.14 mgal. 

Land station navigation was less complex and far more 
accurate. Six of the ten stations were selected for occupation 
because their exact location was clearly shown on the charts 
used. While the other four station locations were transferred 
Moe the chart from hand drawn maps, it is certain that their 
position accuracy was at least twice as good as that of the 


mecanitc stations. 
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2. Measurements 

During the survey the meter vas first lowered to the 
sea surface and a depth counter value vas obtained by nulling 
the appropriate galvanometer; this was repeated as soon as the 
meter reached the sea floor. Next the leveling mode was 
selected and required flood and tilt checks were made. (Only 
OE the 82 bottom stations did the greater-than-15% bottom 
slope indication appear, requiring repositioning of the meter.) 
After the high-speed leveling switch was activated to make 
coarse leveling adjustments, the meter was put into the read 
mode. Tn this mode, while fine leveling adjustments continued 
automatically, the beam position and gravity counter switches 
were alternately used to stop the motion of the beam position 
galvanometer needle in the null position. During this process 
the correct counter value was always approached from the same 
Gmrection (namely, from too low a counter value) to avoid 
hysteresis errors. When vhat was determined to p thnemeeorrecte 
counter value was obtained, an over-adjustment of 0.1 counter 
units was introduced to check for a reversal in the direction 
of M sition motion. This method of obtaining the gravity 
counter value is not only much faster than the slope method 
used with the strip chart record of time versus beam speed, 
but also about as reliable (H. B. Parks, LaCoste & Romberg, 


mie personal communication, 1972). After the gravity counter 


Value was recorded the mass was again clamped in the deck mode 


and the meter raised. 
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The meter was normally on the bottom for about 4 or 5 
min, 2 or 3 min of which were spent in the read mode.  In- 
cluding transit time between stations (but excluding transit 
Ec c and from the first and last stations of each track), 
the 82 bottom stations were occupied in a total of 25 hours, 
yielding an average rate of better than 3 stations per hour, 
or about 18 min per station. The deepest stations usually 
required the most time. 

At the shoreline stations the land gravimeter was 
leveled visually with two bubble glasses and manually operated 
Screws attached to the frame. The scale was then illuminated 
with a light svitch and the mass was manually unclamped.  Cor- 
rect counter values were obtained by simultaneously observing 
the beam position through the microscope eyepiece and manually 
adjusting the measuring screw to null it in the reference 
pscution. 

Fortunately, snycronization between the counter read- 
ings of the bottom gravimeter and the control box was maintained 
throughout the A A A Between the August and Septem- 
ber cruises, however, loss of snycronization did occur once. 
This was evident when a routine counter reading at Station 1 
was several milligals from its standard value, and the loss 
required opening the inner gravimeter case to read its counter 
value so that the control box counter value could be reset to 
agree. 

An important set of measurements midway through the 


survey, accomplished by repetitive occupation, established a 
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3.20 mgal difference between Station 1 on the floor of 
Memeerey Harbor and Station E, the USCG pier benchmark. In 
19η to tieing in all bottom gravity readings with a 
reference benchmark value, this result also related values 
recorded by the two different meters with a station common to 
both surveys. 

BE Environmental Effects 

From a climatological viewpoint, the summer months are 
ESSE suited, due to wind conditions, to gravity operations in 
southern Monterey Bay. By observing surface bubble relative 
motion as well as cable angle, the ship's master was usually 
able to maneuver the R/V ACANIA to avoid being blown downwind 
of the desired position directly above the meter. The wind- 
induced sea and swell presented more of a problem to survey 
operations by causing noticeable pitch and roll. This occa- 
sionally resulted in the gravimeter hitting the extended A- 
frame during lowering and raising evolutions. 

Readings taken in shallow depths were often charac- 
terized by a slowly fluctuating beam position in addition to 
the a nal Ἱ σου ος ehe pean Positron Prior to obtaining the 
final reading. This required observation of the rate of 
change of the position of the fluctuation's mean so that the 
gravity counter switch could be used to eventually hold this 


mean in a fixed position. Since the automatic averaging 


capability of the control box was inoperative, the averaging 
process had to be done by eye, thereby introducing an esti- 
mated reading error of 0.1 counter units into the final counter 


values. 
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The fluctuating beam position was, in part, the visual 
Eucution of a form of micro-seismic bottom motion caused by 
surface swell and wave activity, its influence on beam position 
being a function of depth and type of bottom sediment (LaCoste, 
1967). Additionally, a swaying of the cable was probably felt 
by the meter through the rigidly-attached cable “ο 
The effect of this problem on survey operations was the need 
for more time in the read mode to obtain the correct counter 
value, and, hence, the occupation of fewer stations per hour. 

Fog conditions often prevail in the east central region 
of Monterey Bay during the early morning hours of spring and 
summer. Such was the case during the occupation of Stations 
36, 37, and 38. Its only effect upon the survey was a slight 
HE doration of positioning accuracy; in the absence of 
Visual contact with distant landmarks, radar navigation was 


Ἱ ryzed. 











VO PEA TAPREDUETTON ANDPP RE SENTCADTEON 


The method whereby the raw data collected during the sur- 
vey is transformed into a meaningful CBA picture is discussed 
below. Although much of the theory involved in data reduction 
applies to both land and ocean stations, there are important 
differences. Some of the individual corrections discussed 
AY Only to bottom gravity data; a summary of the reduction 
methods applicable to land stations follows in a separate 
section. 

The majority of the required calculations were accomplished 
BuESUugh use of the NPS IBM-360 computer. Programs were written 
by the author in Fortran language and numerical significance 
was double precision. Copies of these programs are included 


Ferne end of this report. 


ESLLOBSERVED GRAVITY 

As pointed out earlier, gravimeters measure gravity differ- 
ences rather than absolute gravity values. Since gravity 
counter units as recorded from the gravimeter control box are 
numerically close COo but notyexactliy cguiyalent to, mihbligals,; 
they must first be converted to milligals before addition to 
or subtraction from a reference gravity value. The absolute 
reference used for this survey was 979891.7 mgal, the gravity 
value at WH-29 on the Monterey Coast Guard pier. In practice 
a more convenient intermediate reference station at the mooring 
bueys (Station 1) was used to tie in all of the other station 


Beactmgs and tO monitor meter drift. 


43 





The control box counter reading at the pier benchmark was 
3323.05 units; hence, the formula used for computing observed 
gravity at each bottom station was: 

上 29891.7 eve - 3323.05 ) x ( 2.03985 )] (1) 
565 1511 milligals and ος Ce phe Se ieee OD seated) Gravity 
Smee r value at the station in question. The 1.03985 factor 
converting counter unit differences to milligals was obtained 
from a table provided by the manufacturer listing (a) a milli- 
gal value associated with each even hundred counter units as 
well as (b) conversion factors for values between the counter 
units listed. Since all of the stations' observed counter 
values fell between 3300 and 3400, only one such conversion 


τος from the listing was required. 


Pee GALTLTUDE CORRECTION 

Mie ti1rst correction (although the order of various cor- 
rections is not defined in any mandatory way) applied to the 
observed value of gravity 1S normally the subtraction of the 
theoretical value of gravity based on geographical position. 
The determination of theoretical gravity obviously requires a 
M eise ea@owledge Of the shape of Ene Carin so that the dis— 
mance from the surface to the earth's center and the centrif- 
ugal force are known for any location. The actual geometry 
has been determined to be very nearly a triaxial spheroid, or 


ellipsoid of revolution having depressions along the two 45° 


Etude lines and a flattening and bulging of the equator 
(Dobrin, 1960). Years of geodetic investigation have yielded 


numerous revised coefficients for the equation of this 
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reference spheriod to more accurately approximate the earth 
geoid. The author's selection of a particular equation was 
governed (as were all aspects of data reduction) by a desire 
Eemeontorm with the procedures of USGS and the California 
State Division of Mines and Geology in an effort to render 

the final CBA map compatible with existing land gravity maps. 
The coefficients used were those of the 1930 International 
Eerord, the theoretical value of gravity being a function of 
ExErude only (Dobrin, 1960): 

Ją = 978049.0 (1 + 0.0052884 sin?L — 0.0000059 sin?2L) (2) 
aere L is the latitude. Equation (2) results in a predicted 
Freuen to-north increase in the value of local theoretical 
ΙΙΙ Of 1.44 mgal/n mile. 

The assumption that the geoid and the reference spheroid 
coincide would appear erroneous at first glance, since, in 
EU they can differ vertically by as much as 50 m in some 
RES Ons of the world. However, a more meaningful measure of 
ers walidity of ee of this assumption is the change 
in this difference over the area of interest. In the small 
domain of Monterey Bay, the change in the distance from the 
equipotential surface of the geoid to that of the reference 
spheroid is on the order of a few centimeters. It may there- 
fore be concluded that final CBA values will still result in 
an accurate representation of small-to-intermediate scale 


geological features of the upper crust (Grant and West, 1965). 
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MARTH TIDE CORRECTION 

EUSthat the interior of the earth possesses finite 
elasticity, the attractive forces of heavenly bodies (partic- 
ΠΥ the moon and sun) continually act to deform its shape. 
As is the case with the ocean surface, the pertinent manifes- 
tatron of this phenomenon is a small scale vertical fluctuation 
Seene earth's crust which can affect observed gravity values 
to a measurable extent. The true distance from a point on the 
aCe to the earth's center can change by as much as 1 ft 
in a matter of hours, representing a 0.1 mgal change in gravity; 
thus, the effect of these earth tides on gravity measurements 
made over such an interval of time must be considered. 

mie USGS earth tide Fortran program was utilized to obtain 
the required corrections. Dates and times of station occu- 
Mean along with a geographically-central location of 36? 
a2" N, 1219 52' W were the input parameters. The output values 
were tabulated at 20 rin intervals. A copy of this program may 


be obtained from the NPS Department of Oceanography, Code 58Ad. 


Bee DRIFT CORRECTION 

Feadings taken at Station 1 before and after each track 
meme, more often than not, slightly different. These differ- 
ences were attributable to three causes: 


剛三 ar も tidal variation during the time interval 
between base station occupations; 


(5) ocean tidal variation during the same period; 


E meter drift. 
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Since the ACANIA was secured to fixed mooring buoys in an 
east-west orientation, possible north-south variation in the 
mal location of Station 1 before and after each track was 
MmegeeGible. To determine true meter drift it is only neces- 
Sary to remove the effects of (a) and (b) above from the total 
"apparent drift" actually measured. Recalling that counter 
unit differences are very nearly equivalent to milligal dif- 
ferences, earth tide variation was removed directly from counter 
EUungs. Similarly, the changes in attraction due to differ- 
ent sea surface tidal heights above the meter were removed. 
(Each 1 ft of water above the meter decreased observed gravity 
2020727 mgal.) The resultant drift values for each track 
were then applied linearly over the range of stations. The 


greatest drift rate was 0.036 mcal/hour. 


Pee RAE —ATR CORRECTION 
The formula for the value of the Newtonian portion of 


EEUU (in milligals) for a unit mass at sea level (gp) 1s: 


= 


== wo 
gn - R2 

where G is the universal gravitational constant (6.670 x 
10”cm’g”lsec” 2), M is the mass of the earth, and Ris the 
Etance to the earth's center. The free-air correction 

ENG) to gravity data accounts for the fact that the measure- 
ουδ was made at some elevation other than mean sea level, the 
assumed height of the surface of the reference spheroid. (A 


gravity station located exactly at mean sea level, then, would 


Mot require this correction.) Disregarding the mass between 
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1 


the station and sea level for the time being, the free-air 
Eorrection is obtained by differentiation of (3) (Heiskanen 
and Vening Meinesz, 1958): 


ddn 


== Ec Xx t 
FAC Hs (H) 


«715 the elevation of the station above mean sea level. 


(0.09406 mgal/ft) x (H) (4) 


Aesimple modification of this equation yields a formula which 
Eres to bottom gravimetry: 


FAC = (0.09406 mgal/ft) x (Water depth - Tide height) 
(5) 


Since equation (3) implies decreasing values with increasing 
w ia の GTTOm Ehe center of the earth, the free-air correction 
is negative for underwater stations. 

ENigrst-order gravity anomaly, the free-air anamoly (FAA), 
EN  -3:ned with the inclusion of this correction along with 
others mentioned thus far: 

lp oe + BL + D = FAC (6) 
MSIE do, θε, and FAC are as defined by equations (1), (2), 
es), respectively, ET being the earth tide correction, 
BED the drift correction. 

The free-air correction effectively refers station level 
to sea level, and is easily the largest single correction to 
mewapplied. It is therefore clear from equation (5) that 
accurate values of station depth and tide level are critically 
important if an acceptable degree of precision is to be 
achieved. The pressure sensor in the gravimeter chamber pro- 
vides depth counter indications from the sea surface to the 


sea bottom; the difference in the surface and bottom readings 
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converted to depth in feet by means of a linear relation 
extracted from sensor calibration graphs provided by LaCoste 
and Romberg. 

All tide heights used were based on the assumption that 
the existing tide level at each station in the survey area 
was identical with that simultaneously recorded in Monterey 
Harbor. Levels were recorded by the NPS tide recording de- 
wee located in the assistant harbormaster's office on 
Municipal Wharf No. 2 in Monterey Harbor. A reference line 
was drawn on these records at 5.5 ft above staff zero. Mean 
sea level is 5.9 ft above Monterey's staff zero (H. V. Maixner, 
personal communication, 1973). Thus, tide height distances 
were measured relative to 0.4 ft above the reference line. 
These tide values were compared with appropriate values from 
ethe published 1972 tide prediction tables (U. S. Department of 
Commerce, 1972); agreement was very good. In order to make 
mies comparison it was necessary to determine the vertical 
distance from mean lower low water to mean sea level for 
Monterey, since all tide table values are referenced to mean 


lower low water. "ume dilference wac ound tosbo 2.0 NEÉt. 


OR BOUGUER CORRECTION 

While the free-air correction is concerned with changes in 
the denominator of the right-hand side of equation (3), the 
Bouguer correction (BC) addresses changes in the numerator. 
This correction considers how the attraction between the mass 


within the gravimeter and the earth's center of mass is 
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affected by (a) the presence of material between the station's 
elevation (or depth) and sea level, and (b) the density of 
this material. 

The method of the Bouguer correction first assumes that 
the measurement was made on an infinitely flat ocean bottom, 
so local eun irregularities are temporarily neglected. 
Removing tne upward attraction of an infinite slab of over- 
lying water vhose thickness equals the station depth would 
effectively replace that slab ον ασ aursondsresulerin 
Ane: value of observed gravity. Since the formula for the 
ESI on of such an infinite slab is 21TTGoh, where Q is 
BH Eensityv and h is the thickness (Dobrin, 1960), the first 
Passo: the Bouguer correction (BCı), the removal of the sea 
water's ως Der 

BC, - 2TG$, Z (7) 
where Z is the pressure sensor depth and q, is the density of 
sea water. The second part of the Bouguer correction (BC9) 
τ for the increased downward attraction for the sea 
water that would have existed had the measurement been made at 
mean E EV (thus the previous inclusion of the free-air 
correction is assumed): 

BCó - 2]TG Oy (2-24) (8) 
EE 1. Zis the tide height relative to mean sea level. The 
Bard and final part of the Bouguer correction (BC3) considers 
the additional downward attraction that would have existed had 
the material between sea level and the ocean bottom been an 


Mes lab Of crustal material instead of water: 
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ο uoo o) (Z-2.) En 
ESL. Ss the density of average crustal material. Figure 12 
Summarizes these considerations schematically. The total 
Beua@uer correction is then: 

BC = 211 GQyZ2 + 21TG Q, (2-24) + 21PG (Pr -Pw) (2-2) (10) 
SACO mbining terms, 


BC = BC) + BCy + BC3z = 2TTGp,2 + 2 TG op, (2-2 (11) 


t) 
The density of Sea water was taken as 1.027 g/cm?, and an aver- 
age crustal density of 2.670 g/cm? το" ο tor or. Tt is Clear 
maomerigure 12 that the Bouguer correction is positive for under- 
E Stations. Once it is included with the other corrections, 
the simple Bouguer anomaly (SBA) emerges: 
SE Se + ET + D- FAC + BC (12) 
E TOM equation (6), 
SBA — FAA + BC. (ES 
Just as bottom gravity data is reduced to the complete 
Pemguer anomaly in order to compare it with land gravity 
anomalies, its reduction to a gravity anomaly intermediate 
Eg cen the free-air and simple Bouguer anomalies will permit 
comparison with FAA values as determined by a Surface ship. By 
emily adding Ber dabo m om eae ar anomaly of the bottom 
data, such a value ds een ο ατα "τες. his may 
be referred to as a "mass-adjusted" free-air anomaly (FAA'). 


amne he frece air and BonBguer corrections are both func- 


tions of station depth, pressure sensor errors must be exam- 


Meco determine their effect on these corrections. A l-ft 


BHe]puMerror produces a combined FAC/BC error of 0.05 mgal 
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Figure 12. Schematic Representation of the Free-Air 
and ger Corrections: (All densities are 
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(the two corrections are of opposite sign). Although the 
manufacturer's claim of pressure sensor accuracy is 0.5$ 
Es would correspond to + 2.3 ft for the deepest station), 
a more conservative assessment was used. It was assumed that 
ρε} determination was in error by 4 ft at most, so the 


combined FAC/BC error is t 0.20 mgal or less. 


fee fe RRATIN CORRECTION 

Easily the most time-consuming aspect of gravity data re- 
Burn is the topographic, or terrain correction (TC). Often- 
times this correction may be ignored altogether (daD: 
Rietman, personal communication, 1972), and this has been 
generally true for previous underwater gravity work (much of 
which seems to have been done in the Gulf of Mexico). However, 
due to the close proximity to Monterey Bay of very deep abyssal 
plains in the Pacific, as well as such intermediate-scale 
features as the Monterey Submarine Canyon and various coastal 
Mountain ranges, this sort of simplification would be unreal- 
uc at best.  (In.the final result, quite noticeable differ- 
Pees were found in the positions of SBA and CBA isolines 
Spaced at an interval of 5 mgal.) Consider a gravity station 
on the ocean floor located on the upper rim of a deep trench. 
The fact that the trench is filled with sea water produces a 
smaller value of observed gravity than would be measured if 
the trench were filled in with solid crustal material. The 
upward attraction of a nearby guyot would also Lae to decrease 
observed gravity. Thus any terrain deviation above or below 


Station level reduces observed gravity relative to its flat- 


53 


bottom value, both on land and under the sea surface. This 
correction is therefore always positive. 

The generally accepted method for calculating the 
gravitational attraction of topographic irregularities was 
first developed by Hayford and Bowie in 1912. It consists of 
EXE mating the volume of excess or deficit terrain surround- 
mag the station with a series of concentric cylindrical shells 
of varying height (depending on the local elevation or depth). 
eplan view this appears as a set of circular zones; the 
zones are lettered alphabetically from the station outward and 
divided into many compartments. This bulls-eye representation 
is then put on glass or acetate templates scaled to conform 
fee appropriate topographic or bathymetric charts. By 
centering such a template on the station, the average elevation 
or depth of each compartment is estimated to determine its dis- 
tance above or below station level. This difference between 
compartment and station level is the entering argument for 
tables which have been developed to give the corresponding 
EESGV)tational attraction in mEtligals, Mine scompartment cor: 
rections (199 in all) are then summed to yield the total 
terrain correction. 

The tables used in this research are modifications by the 
USGS of the work done by Swick in 1942. His tables are based, 
ENCUurn, on Bullard's 1936 modification of the original Hayford- 
Bowie paper that presented the combined effects of topography 
and isostatic compensation. These tables assume a removal/fill- 


in density of 2.670 g/cm?, and give a 0.615 multiplication 
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factor for the corrections in oceanic compartments. (Since 
EE but sea water of density 1.027 g/cm? fills oceanic 
Gemeartments to sea level, only 1.643 g/ cm? additional mass 
BE -utv, or 61.5$ of 2.0670 g/ cm? , is needed to complete the 
errain Correction density fill-in.) As one would expect, 
the tables' corrections decrease for a given elevation dif- 
ference with increasing distance from the station, and increase 
for a given zone with increasing elevation differences. Thus 
EL -vatron difference effectively adjusts the value of the 
numerator of equation (3), while the zone specification does 
the same to the denominator. 

The effect of terrain upon a gravity value measured on 
the sea floor is more complex than its terrestrial counterpart; 
Es understandable in view of the inclusion of an additional 
el the sea water, surrounding the station. This effect 
must be computed in three stages; Figure 13 presents the analy- 
Sis in cross-section. 

bey the Sea floor (Fig. 13-A) been horizontal, a terrain 
correction would ο be required. The existence of a uniform 
ENEGD of 1.027 g/ cm? material extending infinitely in all 
directions ech station level and sea level can be neglected 
without changing the TC value to be calculated since the topo- 
graphic effect of such a slab of given density and of constant 
thickness equals zero. The first step shown in Figure 13 
simplifies the problem and requires no actual calculations; 
it effectively replaces the sea water in Sector 1 with air and 


decreases the crustal density in Sector 4 to 1.643 g/cmJ. The 
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Etant mass distribution surrounding the station shown in 
Eure 13-B has an attractional effect (due to terrain) iden- 
σα to that shown in Figure 13-C, i.e., as if the station had 
been located right at the coastline of a layered continent of 
two different densities. 

Peco two Consists of a “normal” terrain correction as de- 
Scribed earlier. Land compartment corrections are sub-totaled 
added to 61.5% of the sub-total of the corrections for 
oceanic compartments. This step effectively transforms the 
sea water in Sector 2 into rock of density 2.670 g/cm, and 
removes rock of similar density from Sectors 4 and 5. Sector 
5 1s now air (as desired), but the use of a density of 2.670 


3 was excessive in Sector 4, since its density was only 


g/cm 
1.643 g/cm3 νετ step l (Fig. 13 D): 

Eco three corrects this situation by "putting back" 
material of density 1.027 g/cm?, thereby making Sector 4 a 
Boon of air. This last step is accomplished by subtracting 
from the previously accumulated total correction after step 2 
a correction based on only one elevation (station depth) fer 
all ara compartments. It is 38.5% (the density ratio of 
EIN ter to crustal rock) of the subtotal thus obtained from 
the tables. 

Adding the total terrain corrections to the data refers 
the resulting gravity anomalies to their values had the station 


Been located on a flat earth crust (Fig. 13-F). Table II 


summarizes the steps just described. 


a 












Seclor 





Designation 


Station 









Step Three 


ἭΝ 





Decrease 



















OZ 
g/cm? to yield 
zero effective 
density. 







present by 


an amount 
2.670 g/cm3. 


equal to sea 
water density 















3 
No Change No Change. 


Table II. Summary of Terrain Correction Procedure for a 
Gravity Station on the Sea Floor. (Steps one 
and three are unique to bottom gravimetry.) 






Decrease 
rock density 
EO Zero. 


1 sea water 
density to 
zero. 
Increase 
M ge T 
Sector 2 No Change. density No Change 
τν ΕΠ RO 
EOCK. 
= A u μμ 
Decrease rock| Decrease m 
; density by density Present, densı- 
5 






58 


= - 





It would not be possible after step l to compute separate 
EsmEections for Sectors 4 and 5 with density ratio factors of 
MES ana 1.0, respectively. This is due to the: fact that 
the tables were developed with the assumption that the base of 
each land compartment is at station level. 

Most of the maps, charts, and templates required were 
provided by USGS. One such template was constructed by the 
author, however, scaled for use with U. S. Coast and Geodetic 
ο ου (USC&GS) Chart No. 5403. The bathymetric information 
on this chart was considerably better than that of others 
available, and was especially helpful for corrections in Zones 
τουσ I. 

ime Should be kept in mind that the slope of the bottom in 
Figure 13 is greatly exaggerated. Since bottom slopes in most 
of southern Monterey Bay are on the order of 1:100, many 
Stations' cumulative corrections as far out as Zone G totaled 
zero. For this same reason the considerations of steps 1 and 
NEL n practice, not required. Sector 4 compartment ele- 
vations could never exceed station depth and were always 
negligibly smaller within the closest several zones, with the 
exception of Station 65 (located on the eastern end of the 
thalweg of the Monterey Submarine Canyon). Here, a simple 
we pl ication by the 0.615 density ratio factor sufficed for 


the Sector 4 corrections. As this number is usually used to 


Enssorm water in Sector 2 into rock, it similarly transforms 
rock in Sector 4 into water. This simplification was possible 


because Sector 5 corrections for Station 65 are nill (since 
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the elevations of the area surrounding Moss Landing are so 
small). The fact that the Sierra Nevadas are just beyond 
the 100-mile radius of the most distant zone (Zone 0) and 
that the San Juaquin Valley is so UST TON Plat wend close to 
sea level also served to expedite calculations somewhat. 

A realistic assumption regarding possible errors in 
Eating compartment elevations is + 0.02 mgal per zone. 
If errors of this magnitude and similar sign existed in each 
Berne 15 zones, the total terrain correction error would be 
EE 50 mgal for each station. 

mie rOregoing analysis treats the terrain correction in- 
dependently of the earlier Bouguer correction. This is valid 
Since the Bouguer correction assumes a flat station environ- 
ment (Heiskanen and Vening Meinesz, 1958). Independent analv- 
Sis is further justified by considering a change in the order 
@eemMaking the corrections: applying a pure terrain correction 
Bothe free-air anomaly, followed by the Bouguer correction, 
ExUdoresult in an identical CBA. Regardless of order, the 
free-air correction only looks at vertical distance, the 
Euer correction deals with vertical mass attractions, and 
the terrain correction accounts for quasi-horizontal mass 


Ferractıons. 


BEZ CURVATURE CORRECTION 
The curvature correction (CC) may be thought of as a 


terrain correction of macroscale proportions since it accounts 


τ tne curvature of the earth. (Recall that the assumed 
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Mguer slab was infinitely flat.) For shallow bottom stations, 
ES correction is negative and is given by equation (14) 
mu. Robbins, personal communication, 1972): 


ιτ. 1422 x 107% H - 3.282 x 1078 H2 + 1.270 x 1015 3 
(14) 


lere CC is in milligals and H is height in feet above sea 
vel. (Oceanic depths are entered as negative numbers.) 
When the terrain and curvature corrections are applied to 


le simple Bouguer anomaly, the result is the complete 


(|uiguer anomaly (CBA). Its value may be stated in equation 
Nm: 
EN «9. * ET * D- FAC * BC * TC - CC (15) 


E from Ίου 2): 

CBA = SBA + TC - CC 22016) 
le CBA is the difference between observed gravity and the 
ique expected had the station been located on a flat solid 
irface at sea level, unaffected by astronomical attractions. 
i18 final gravity anomaly refines the geological implications 
- tne SBA and represents a useful tool that geologists can 


EE ρο to infer substructure mass distribution. 


SHORELINE SURVEY < 

Data reduction for the ten stations ashore proceeded in a 
Eon Similar to that of the bottom survey, but with five 
ceptions. First, since calibration factors for the con- 
»rsion of counter units to milligals are unique to each 
irticular gravimeter, the equation for Jo for the Model G-17B 


-avimeter was not the same as equation (1). Second, free-air 
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corrections were added (rather than subtracted), since on 
land, mean sea level in general is closer to the earth's 
enter than station elevation. Third, Bouguer corrections 
were subtracted (rather than added), because part of the 
observed attraction was due to the material between station 
level and sea level. Thus, while underwater Bouguer correc- 
tions are composed of three parts, terrestrial Bouguer correc- 
tions need compensate for only a single attraction. Fourth, 
the three-step terrain correction analysis does not apply. 
Finally, since all elevations were within 10 ft of mean sea 
level, curvature corrections were not required. (A station 
elevation of 22 ft above or below sea level is required to 
pa duce a curvature correction of 0.01 mgal.) The equation 
EEumerizing land station data reduction is then: 

EN G4 7 de + ET + D+ FAC το (17) 
Estimated error values for land station corrections are 


gmoeo2 mgal for g tO 0d MG 2 σι ο cH C pon be 


O 1 
Epssutroning accurate to within 1/20 n mile), + 0.06 mgal/ft 
heeestation elevation, and, again, + 0.30 mgal for the terrain 


memerection elevation estimations. 


J. DATA PRESENTATION 

The values of the "mass-adjusteda" free-air anomaly and 
complete Bouguer anomaly for both shoreline and bottom stations 
were recorded on charts next to appropriate station locations 
and scalar analyses were performed by hand with isolines 
spaced at a 5-mgal interval (Fig. 14 and 15). Table III 


presents the values of these and two other gravity anomalies 
62 


A — ο ος 0 — 








of interest for the 92 stations; values of observed and 
theoretical gravity as well as values of the various correc- 


tions are presented in Appendix A. 
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Figure 14. "Mass-Adjusted" Free-Air Anomaly Map of Southern 
Monterey Bay. (Contour values are in milligals.) 
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Figure 15. Complete Bouguer Anomaly Map of Southern Monterey 
Bay. (Contour values are in milligals.) 
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VI. DATA ANALYSIS AND GEOLOGICAL INTERPRETATION 


Figure 15 is in general overall agreement with the 
regional trend of CBA values on adjacent land areas as deter- 
mined by Fairborn (1963), Sieck (1964), Bishop and Chapman 
NEC. and Ivey (1969). The shape of the +10 and +15 mgal 
isolines is in excellent accord with the shape of the granitic 
εξ north and east of Pt. Pinos (Fig. 3). Additionally, 
there is a discernible gradient increase from the +10 to -5 
mgal isolines in the southernmost region of the Bay; this 
seme iaes with the Tularcitos Fault Zone. The previously 
assumed existence of -30 mgal isolines seaward of the Marina 
area (Chapman and Bishop, 1967) was due to insufficient station 
coverage along the coast and no offshore coverage; the -25 
mgal isoline in Figure 15 is confirmed by a CBA value of 
-25.620 at Station H. Monterey Peninsula CBA values were 
mma tO be slightly higher than those presented on the 1967 
Santa Cruz Sheet, but, again, bottom gravimetry CBA values 
were E fied by the independent land survey (Stations A, B, 
C, and D). Of far greater geological significance than the 
actual CBA value associated with each isoline is the general 
shape and trend of these lines, an aspect of Figure 15 that 


all investigators concerned could easily accept. 


Possible error magnitudes are summarized in Table IV. 
These estimates are considered to be maximum values; it is 


felt that the final CBA values of the author's survey are 


Beeurate to within 1 mgal. 
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ENDS readily apparent from Figure 15 that the north- 
eastern portion of the survey area is characterized by a 
deficiency of mass density relative to the southwestern 
region. This deficiency represents an increase in depth to 
the top of the dense granitic basement from southwest to north- 
east, and is in very good agreement with seismic reflection 
profiles presented by Greene in 1970 (Fig. 16 and 17). Since 
a single geophysical method seldom assures a unique solution, 
this confirmation of previous seismic reflection interpretation 
Ex numportant result of the author's study. The CBA pro- 
uu in the figures lack the detail (particularly in the 
Tularcitos Fault Zone) of the corresponding seismic profiles; 
meets most probably due to the grid size of the station net- 
Lo Additionally, there may exist an insufficient density 
Sem@erast, although this possibility is remote in view of the 
apparent displacement and the density difference between the 
Monterey Formation and granodiorite below Monterey Bay. 

Figures 16 and 17 indicate an average complete Bouguer 
gravity anomaly decrease of 5 mgal for each 520 ft increase in 
depth ES the granitic basement structure. Since the 11 n mile 
extent of profile A-B and corresponding 3700 ft drop in the 
level of the basement represent a slope at the interface be- 
tween the Monterey Formation and the granodiorite of only 


about 3?, we can approximate the difference in attractive 


forces at two points along that profile whose CBA values differ 
by 5 mgal by the attraction of a horizontally infinite slab 


520 ft thick whose mass density is equal to the density 
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Figure 16. Comparison of CBA and Depth of Granite Substructure 
as Determined by Seismic Reílection for Profile A-B 
(Upper profile after Greene, 1970). The ¡location 
of Profile A-B is shown in Figure 15. 
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Escure 17. 


Comparison of CBA and Depth of Granite Sub- 

structure as Determined by Sewsmic Reflection 
for Prortilegcp (UmsbessoromMcudhbcer Greene,u1970). 
The location of Profile C-D is shown in Figure 15. 
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difference between the two rock types. So doing permits 
ENuating a theoretical attraction of 2T G p. h (where h is 
Beet and pq. is the density contrast) with the associated 

5 mgal inferred value, and yields a density contrast of 0.75 
g/cm. πο the slope of the igneous/sedimentary rock inter- 
face is fairly constant along the profile, the values of sub- 
layer terrain irregularity corrections at two such locations 
will be close and their difference can therefore be neglected. 
This analysis assumes, then, that the relation between CBA and 
depth to basement is due solely to the magnitude of the 
density deficiency existing in the upper formation; the 
attraction at the point farther toward the northeast will be 
less since a thicker layer of the (lighter) Monterey Formation 
exists below. 

A fairly accepted value for the density of the granodiorite 
m. 75 g/cm?, but the density of the in situ Monterey Forma- 
tion is not accurately known. While this formation consists 
of multiple members, the Miocene marine shale predoninates; 
its B suci aensity was determined experimentallv by 
the author to be 1.85 g/cm? (true in situ density is higher 
Since the lab samples were not under pressure).  Sieck (1964) 
found the dry density of Monterey shale to vary from 1.41 to 
220 g/cm?, with an average of 1.80 g/cm?, and his measured 
density of Monterey sandstone averaged 2.10 g/cm>. Exami- 
nation of various references to shale densities in general 
indicates that it is a highly variable figure (depending on 


age, location, and many other factors); values for diatomaceous 


7 3 








shales in the San Jauquin Valley are between 0.9 and 1.1 g/cm, 
while some water saturated shale types can be as dense as 

3.21 g/cm? (Jakosky, 1950). In view of these many factors, 

the 2.00 g/ cm? implied density of the Monterey Formation (2.75 
g/cm? Granodiorite minus the gravity survey's calculated 0.75 
g/cm? density contrast) would appear quite acceptable as a 
mue” value. 

Having arrived at densities of 2.00 and 2.75 g/cm3 for 
the two major geological units in the area, it is now appropri- 
ate to consider the validity of using 2.67 g/ cm? for the 
assumed average crustal Bouguer density. Since the Monterey 
Formation is less than 1 km thick, while the granodiorite sub- 
EE cture is probably several times thicker, the average 
ου αι density in the vicinity of southern Monterey Bay is 
IE loser to 2.75 than to 2.00 g/cm?. Thus, while the 
frequently used average crustal density of 2.67 g/ cm? can be 
unrealistic in many regions of the world, it is probably not 
the case here. 

Throughout the course of this research it was anticipated 
that Tine positions over the deep Monterey Submarine Canyon 
ENS be inferred from the results of the bottom surveys in' 
the northern and southern halves of Monterey Bay. Both sets 
of data were obtained with the same equipment in an identical 
manner and reduced in a comparable fashion, so there is no 
meason to doubt the validity of making such an inference. 

ο. το 18 is the result of combining the data from Cronyn's 


(1973) survey with that of the author. 
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Figure 18. Composite Gravity Map of Monterey Bay (Northern 
Monterey Bay contours after Cronyn, 1973). 
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The fact that the southernmost positions of Cronyn's knovn 
E50 -10, and -15 mgal isolines must bend to meet the northern- 
most positions of those of the author supports the theory that 
a fault exists along the axis of the Monterey Submarine Canyon. 
The difference in lithologies on opposite sides of the Canyon 
led Martin and Emery (1967) to recognize the so-called Monterey 
Fault parallel with the Canyon and offset at numerous locations 
by transverse faults of different magnitudes. Although the 
fault's suggested orientation is at an angle to the trend of 
almost all of the many fault zones in central California, the 
proposed situation is seen elsewhere, as in the case of the 
Santa Ynez Fault region north of Los Angeles, where the strike 
of the San Andreas Fault changes. Further evidence of the 
existence of the Monterey Fault based on seismic reflection 
tracks across the Canyon has been presented by Greene (1970). 

Although some of the advocates of the Monterey Fault 
suggest that the motion along the fault is dip-slip with the 
northern side having descended relative to the southern side, 
this is not uniquely indicated by the gravity data. On both 
Fides of puescanyone (rin the eastern on mOr themBay) the 
CBA isolines run north-south with the higher values seaward 
(Fig. 18 and 19-A). If the motion of the Monterey Fault is, 
Mmract, primarily dip-slip, the™shape” of the CBA field over 
the Canyon would suggest a relative motion such that the 
southern side has descended rather than the northern side, as 
πη in Figure 19-B. This is supported by the fact that the 


motion of the northeastern side of the Tularcitos Fault is 
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Figure 19. Two Possible Structural Explanations for the 
Inferred CBA Field in the Vicinity of the 
Monterey Submarine Canyon. 
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E uard (Greene, 1970). If, on the other hand, the fault 
muon is primarily right-lateral strike-slip, the inferred 
CBA isolines over the Canyon may be explained as in Figure 
19-C, wherein the northern side has moved east relative to the 
Southern side. Specifying which of the two types of relative 
pon has occurred is not possible from the gravity infor- 
mation alone; quite possibly a combination of both types of 
Me~eron has occurred, causing an oblique fault. 

It may well be that the intersection of the Monterey Fault 
ememehe Palo Colorado - San Gregorio system farther to the 
west corresponds to the area where a concentration of recent 
seismic epicenters is located. More geophySical and ocean- 
ographic study of this region is needed before an accurate 


Structural model can be developed. 
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Vee ore WORK 


te is recommended that CA-259 (the Monterey County Airport 
station) of the state gravity base station network be re- 
occupied in accordance with procedures outlined by Chapman 
(1966) to re-establish a correct value. Local construction 
there will be completed soon. 

[emis suggested that additional gravity measurements be 
Made in the vicinity of the Monterey Submarine Canyon. Where 
a large portion of the coverage needed might be precluded 
(with a bottom gravimeter) by excessive depths and/or slopes, 
Surface ship gravimetry could be called upon to complete this 
task. One such survey consisting of several canyon crossings 
was conducted by USGS in November, 1972, aboard the USNS 
BARTLETT, but at this writing the results are not yet available 
for inclusion here. They will eventually help greatly in re- 
moving the question marks in Figure 18. 

Bthis writing the gravimeter and auxiliary equipment 
used for the Monterey Bay surveys is back aboard the R/V 
ΙΙΙ. The author has thoroughly familiarized two new NPS 
students with the entire operation, and surveys are underway 
to extend shallow water coverage in two areas, one from Pt. 
Bonos to Pt. Sur and the other from Santa Cruz northwest to 
beyond Davenport. The goal of the project sponsors is to tie 
in these surveys with those conducted by the author and co- 
workers, so that a continuous picture of the off-shore Bouguer 


anomaly will become known. Hopefully this goal will be achieved. 
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APPENDIX A (Continued) 
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COMPUTER PROGRAM FOR UNDERWATER STATIONS 
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